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We performed point-contact spectroscopy on the binary superconductor MgB2. The differential 
conductance shows gap-related structures which vary in width and position from contact to contact. 
The distribution of energy gaps shows a distinct accumulation around 1.7 and 7meV which is 
associated with the occurrence of a small and a large energy gap in MgB2. While with increasing 
T the structure in dl/dV associated with the small gap is present up to Tc, in magnetic field it is 
suppressed well below Bc2- 
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The observation of superconductivity at 39 K in the 
simple binary compound MgB2 was a surprise for most 
condensed matter scientists and immediately raises the 
question of the mechanism of superconductivity in this 
material. Following a basic approach to superconduc- 
tivity, two key issues have to be considered: the order- 
parameter symmetry and the coupling mechanism. For 
the latter, a significant boron isotope effect was observed 
in MgB2 1^ which is consistent with a phonon-mediated 
BCS mechanism where the boron phonon modes are play- 
ing an important role. A moderatly strong electron- 
phonon coupling constant Aci_ph « 0.8 was derived from 
measurements of the specific heat in good agreement 
with recent theoretical predictions . 

The order-parameter symmetry has already been in- 
vestigated by tunneling spectroscopy and point-contact 
spectroscopy (PCS) fTj-p^. Although the spectra show 
unambiguous features of an energy gap in the density 
of states, probably with s-wave symmetry of the or- 
der parameter, the results are controversal on the gap 
width. Values of 2A/kBTc ranging from 1.2 to 4.7 have 
been reported, raising the possibility of an anisotropic 
energy gap or even multiple gaps. Careful analysis 
of specific-heat measurements |l^,|l^ and photoemission 
spectroscopy Q support this scenario as well. 

We present measurements of the differential conduc- 
tance dl/dV of normal conductor/superconductor point 
contacts between a Ft tip and a MgB2 pellet. The MgB2 
pellet was prepared in an Ar atmosphere from magne- 
sium filings (> 99.9%, Chempure) and powdered boron 
(99%, Ventron). The powders were mixed in the ratio 
Mg:B = 1.1:2 and pressed into a pellet. After wrapping 
the pellet in a Ta foil it was enclosed in an Fe cylinder 
and inserted in a high-pressure furnace. The pellet was 
then heated to 1073 K for 1.5 h and to 953 K for 2h un- 
der an Ar pressure of 58MPa. The DC magnetization 
of the pellet (Fig. |l|) indicates an onset superconducting 
transition temperature r°"'"=* = 38.2 K. Almost the full 
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FIG. 1. DC magnetization of MgB2 in different applied 
magnetic fields. 

theoretical shielding signal is reached in 1 mT applied 
field in a zero-field cooled (zfc) - field heated (fh) tem- 
perature cycle. The Meissner signal remains weak, i. e. 
of the order of a few percent. Upon reducing the applied 
field to 10 /iT, a kink in the transition curve at 35 K re- 
veals the granular structure of the sample: below T""'*^* 
first the individual grains become superconducting, be- 
low 35 K intergranular shielding currents cause shielding 
of the whole pellet. 

PCS spectra were taken on as grown and polished sur- 
faces of the MgB2 pellet. Different contacts were es- 
tablished in situ at low temperatures in liquid Helium. 
The PCS spectra, i.e. the differential resistance dV/dl 
vs V were obtained using a standard four-point lock-in 
technique. The zero-bias resistance Rq of stable contacts 
ranges from 1 to 60 fi. dl/dV vs V curves are obtained 
by numerical inversion of the dV/dl data. The dl/dV 



1 



10 - 



i 5h 







-i — I — I — I — I — I — I — I — I — I — I — I — I — I — |- 



MgBj - Pt 



BCS value 



. nn n . 



EL 



0123456789 
A (meV) 

FIG. 2. Extracted values of the gap energy for different 
MgB2 contacts 

curves taken at low T <^ (between 1.7 and 5K, cf. 
Fig. H), show most often maxima at voltage ±Vb, pre- 
sumably related to the superconducting energy gap A. 
The position of the maxima A = eVb, varies depending 
on the location of the Pt needle on the polycrystalline 
sample where the point contact was established. The 
histogram depicted in Fig. ^ summarizes the position of 
maxima of 20 contacts. In a few cases a single maximum 
centered at y = was observed (Fig.^). For these cases 
the value of Aq = A(T — > 0) was derived from the mod- 
ified BTK fit (see below) and is also included in Fig.||. 
This histogram in nice consistency reflects the gap dis- 
tribution reported in literature. One dominant feature is 
a maximum at about 1.7 meV, a second occurs at about 
7meV. 

Following the BTK theory of transport through 
a metallic microcontact between a normal metal (N) 
and an isotropic superconductor (S), Andreev reflection 
is the dominant transport process at the N-S interface. 
Andreev-reflected charge carriers lead to an excess cur- 
rent which increases the point-contact conductance. De- 
pending on the strength of the interface barrier a char- 
acteristic double-maximum feature is expected in dl /dV , 
with the maxima at ±A/e where A is the energy gap. In- 
elastic scattering in the contact region usually leads to a 
smearing and weakening of the gap-related features. This 
can be empirically taken into account by replacing the 
quasiparticle energy E with iS-fiF, where F parametrizes 
the finite lifetime of the quasiparticles . We observed 
a few contacts where the dl/dV curves show maxima at 
±Vb and ±V^i. These curves can be regarded as a super- 
position of two double-maximum features. Consequently, 
the maxima are interpreted as being associated with two 
different energy gaps and the presence of multigap super- 




FIG. 3. Measured dl /dV curves and BTK fits to the small 
gap structures for different MgB2/Pt contacts at T = 1.7 K 
(panel a-c) and 5 K (panel d and e) either with thermal smear- 
ing only (r = 0, solid lines) or with additional finite-lifetime 
smearing (F 7^ 0, dashed lines). The fit parameters are sum- 
marized in Table| 



conductivity in MgB2. 

Multigap superconductivity in MgB2 has recently 
theoretically been suggested by Liu et al. based 
on the electronic structure proposed by Kortus et 
al. |^l|: The Fermi surface consists of two sets of 
sheets, two-dimensional (2D) cylindrical sheets and 
three-dimensional (3D) tubular network. In the clean 
limit, two different superconducting order parameter are 
predicted: a larger one on the 2D part of the Fermi sur- 
face, a smaller one (about 1/3) on the 3D part of the 
Fermi surface. The multi-sheetedness of the Fermi sur- 
face should lead to a substantial difference between the 
in-plane and out-of-plane tunneling spectra. The effect 
of strong defect scattering drives the system into the 
dirty limit where the two energy gaps are mixed and 
an isotropic BCS gap is obtained. Quite remarkable is 
the effect on the temperature dependence of the gap in 
both limits. In the clean limit both energy gaps close at 
Tc while in the dirty limit the resultant isotropic gap is 
substantially suppressed. 

In the following we focus on spectra which show 
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solid lines dashed lines 

T (K) Ao (meV) Z r Ao (meV) T (meV) Z 

(a) 1.7 1.5 0.5 0.368 1.5 0.65 0.57 

(b) 1.7 1.7 0.8 0.079 1.1 1.4 0.9 

(c) 1.7 1.7 0.3 0.305 1.5 0.87 0.47 

(d) 5 2.1 0.615 0.128 1.3 1.3 0.73 

(e) 5 - - - 1.1 3.6 0.1 

TABLE I. Fit parameter of the theoretical curves in Fig.lsl 



features associated with the small gap structure only. 
Fig. ||a-e show a number of spectra of this type together 
with theoretical curves calculated within the modified 
BTK theory either with thermal smearing only (F = 0, 
solid lines) or with additional finite-lifetime smearing 
(F 7^ 0, dashed lines). The fit parameters of these 
curves are summarized in Table |. Although the width 
of the spectra is almost the same, the voltage depen- 
dence is different. Curve a and d (referred to as type 
1) are described very well with a pure (F = 0) BTK 
fit with an energy gap (assumed to be isotropic) Ao = 
1.5 meV (a) and Aq = 2.1 meV (d), a barrier strength 
Z ^ 0.5 (a) and Z = 0.615 (d), respectively. The 
fits have been scaled by a factor of r — [(i//(iF(0) — 
d//dy(15mV)]oxp/[d//dF(0) - d//dF(15 mV)]thco = 
0.368 (a) and 0.128(d), respectively, in order to meet 
the absolute magnitude of the experimental conductance 
change [d//dy(0) — c?//(iF(15 mV)]cxp- The common fea- 
ture of curve b, c, and e (referred to as type 2) is that 
the decrease of the conductance for V > A/e is weaker 
than expected from a pure BTK fit, i. e. the experimental 
curves are much wider. Smearing due to inelastic scat- 
tering in the contact region can be one origin of wider 
structures which has been modelled by F 0. The pa- 
rameters A, F, and Z are chosen in such a way that the 
calculated curves describe the experimental conductance 
change without the necessity of scaling. Introducing of r 
as an additional fit parameter to improve the fit would 
require an unphysical value r > 1. A relatively large F 
of the same order of magnitude as A or even larger is 
needed to model the curves, and the resulting descrip- 
tion of the spectra is not really satisfactory. The failure 
of the model in describing curve c hints at a second, in- 
trinsic origin, namely an apparent gap distribution due 
to an anisotropic small gap. This means that in addition 
to the occurrence of two energy gaps in MgB2 , these gaps 
themselves might not be isotropic. Recently Giubileo et 
al. [Til pointed out that the gap width varies even in the 
case where spatial inhomogeneities can be ruled out as 
an origin of a gap distribution. It might be of interest to 
note that all three curves show a weak conductance max- 
imum around zero-bias which might be caused by nodes 
of the gap function jl^ . 

In order to further elucidate the difference of the spec- 
tra we now turn to the temperature dependence of both 
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FIG. 4. Temperature dependence of point contacts of type 
1 (a) and 2 (b). 



types depicted in Fig.^ a and b together with theoreti- 
cal curves. With increasing temperature the gap-related 
structure weakens in both cases. However, while gap- 
related features of the spectra of type 2 vanish close to 
the bulk Tc and almost follow the BCS weak-coupling be- 
havior of A(T), the isotropic gap extracted from spectra 
of type 1 closes at « 0.7 Tc. Without overrating the ex- 
perimental facts, the results resemble the proposed multi- 
gap scenario. Bearing in mind the granular structure of 
the sample and the unknown microscopic structure of the 
point contact, it is conceivable that both clean and dirty 
limits are locally realized on different parts of the sample. 

In an applied magnetic field B the gap structure weak- 
ens with increasing field and vanishes around 6T (see 
Fig.|). The small gap is much more sensitive to an ap- 
plied field than expected from estimates of the upper 
critical field i3c2(0) from resistivity and magnetiza- 
tion measurements |2^]. These data suggest that a field 
well above 12 T is necessary to suppress superconduc- 
tivity at 1.7 K. Szabo et al. have recently reported an 
even stronger field dependence of the small gap |jl^ . In 
line with specific-heat measurments [|l^,|l^ , i? = 1 T sub- 
stantially suppresses the small gap. The result hints at a 
different nature of the two energy gaps probably due to 
their different dimensionality. 
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FIG. 5. Magnetic field dependence of tfie point contact 
sliown in Fig.Hb. Tfie inset sfiows A{B) extracted from the 
BTK fit to tiie data. 

In conclusion, we have examined point contacts on 
MgB2 in order to investigate the superconducting en- 
ergy gap. The gap width obtained from dl/dV vs V 
spectra varies from contact to contact depending on the 
microscopic structure and grain orientation in the point- 
contact region. The result is compatible with multigap 
superconductivity with a small gap present up to Tc- 
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